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ABSTRACT 


The  structure  of  turbulent  wakes  of  hypersonic  bodies  is 
argued  to  be  largely  controlled  by  the  combined  action  of  relatively 
large  scale  random  convective  motions,  which  introduce  inhomo¬ 
geneities  in  the  wake  core,  and  of  molecular  diffusion  aided  by 
turbulent  distortions  of  fluid  elements.  A  previously  developed 
mathematical  model  representing  turbulent  wake  mixing  in  terms  of 
a  mixing  lag  is  discussed  and  qualitatively  justified  in  terms  of  the 
above  described  mixing  mechanism,  and  somewhat  generalized. 
Numerical  predictions  for  the  magnitude  and  variation  of  mass 
density  fluctuations  in  a  relatively  low  speed  wake  are  obtained  by 
solving  a  simpli^’ied  set  of  equations,  in  which  the  lag  is  represented 
by  an  effective  mixing  boundary,  and  are  in  encouraging  agreement 
with  the  general  magnitude  and  trend  observed  experimentally, 
although  firm  experimental  measurements  are  lacking  to-date. 
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1.  INTRODUCTION 


The  study  of  the  structure  of  turbulent  wakes  of  hypersonic 

bodies  has  been  pursued  vigorously  in  recent  years  as  part  of  the 

continuing  research  in  re-entry  phenomenology.  Many  theoretical 

1,2,3 

analyses  of  the  turbulent  wake  ’  *  ’  have  been  based  on  an  extension 
of  laminar  wake  analyses,  through  the  use  of  turbulent  transport 
coefficients  to  represent  the  effects  of  turbulence.  The  proper 
characterization  of  the  variation  of  fluid  properties  and  of  chemical 
reactions  within  the  turbulent  wake  requires,  however,  a  more 
detailed  description  of  turbulent  mixing  than  is  provided  by  analyses 
based  on  eddy  transport  coefficients,  which  are  particularly  defi¬ 
cient  in  representing  the  effects  of  large  local  inhomogeneities  on 
chemical  reactions  in  the  turbulent  core. 

In  an  attempt  to  provide  a  more  physical  description  of  the 

important  features  of  mixing  in  the  turbulent  wake,  a  simple  wake 

mixing  model,  based  on  tb...  known  properties  of  turbulent  convection 

4  5 

and  mixing  in  low  speed  wakes  has  been  previously  formulated  '  . 

In  that  model,  a  distircticn  was  made  between  a  fluid  component 
consisting  of  fluid  ingested  into  the  turbulent  core  and  randomly 
convected  within  it,  but  which  had  not  undergone  appreciable  mixing 
or  diffusion  on  a  molecular  scale  and  a  more  homogeneous  well 
mixed  component.  A  mixing  lag  was  used  to  represent  the  delay  in 
molecular  mixing  of  fluid  newly  ingested  into  the  turbulent  core. 

The  present  paper  reviews  and  discusses  the  physical  basis  and 
validity  of  the  model  in  somev.'hat  greater  detail,  proposes  a  some¬ 
what  more  general  formulation  of  wake  mixing  and  presents  a 
numerical  solution  to  a  simplified  set  of  equations  representing  the 
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mixing  model,  in  which  the  mixing  lag  previously  formulated  is 
represented  by  a  mixing  boundary  which  "lags”  behind  the  turbulent 
boundary. 

The  numerical  solution  was  obtained  for  the  wake  of  a  rela¬ 
tively  low  speed  blunt  body,  representative  of  the  speeds  of  pellets 

fired  in  ballistic  ranges.  The  results  are  in  at  least  qualitative 

.  ,  .  6,  7  ,  ,  . 

agreement  with  experiment  ,  insofar  as  predictions  of  root  mean 

square  gas  density  fluctuations  in  the  wake  are  concerned.  Unfor¬ 
tunately,  the  absolute  magnitude  of  the  fluctuations  observed  experi¬ 
mentally  is  still  somewhat  in  doubt  at  this  time.  Additionally,  it 
must  be  noted  that  while  the  interest  in  hypersonic  turbulent  wakes 
is  centered  on  the  case  of  high  speed  wakes,  mass  density  fluctua¬ 
tion  measurements  have  so  far  been  confined  to  relativ'ely  low 
velocities.  Thus,  the  most  striking  effect  of  the  present  model, 
which  is  its  influence  on  the  wake  chemistry,  and  particularly  the 
ionization  chemistry,  cannot  be  directly  tested  at  this  time.  How¬ 
ever,  the  ability  of  the  model  to  reproduce  the  gross  features  of  the 
mass  density  fluctuations  at  the  lower  speeds  is  an  encouraging  indi¬ 
cation  that  it  may  be  applied  with  some  confidence  to  describe  the 
properties  of  hypersonic  re-entry  wakes. 
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2.  PHENOMENOLOGY  OF  TURBULENT  WAKE 
MIXING  AND  WAKE  STRUCTURE 

The  growth  of  the  turbulent  wake  behind  bodies  in  high 

Reynolds  number  flows,  and  the  processes  of  engulfment  of  the  outt  ■ 

flow  in  the  wake  growth,  and  of  mixing  and  con'  ection  of  the  core 

8,  9 

gases  are  relatively  well  understood  for  subsonic  flows  ’  .  While 

there  undoubtedly  will  exist  differences  between  such  wakes  and  the 
wakes  of  hypersonic  bodies,  it  can  be  reasonably  assumed  that  the 
same  qualitative  features  will  persist  in  the  hypersonic  case,  parti¬ 
cularly  when  it  is  realized  that  the  velocity  fluctuations  in  such  wakes 
are  subsonic  beyond  the  first  few  body  diameters,  due  to  the  very 
rapid  initial  decay  o^  the  mean  wake  velocity.  On  the  assumption 
that  low  speed  turbulence  results  are  applicable  to  hypersonic  wakes, 
the  grow'th  of  tne  barbulent  wake  and  the  manner  in  w'hich  it  engulfs 
outer  fluid  may  be  broadly  characterized  as  follows;  The  local  pro- 
p  gation  of  the  turbulent  front  bounding  the  wako  core  occurs  by 
diffusion  of  vorticity  fluctuations  into  the  outer  fluid  ’  .  The  dif¬ 

fusion  is  initiated  by  the  smallest  eddies,  since  they  repre'^ent  the 
greatest  velocity  gradients.  The  transfer  of  turbulence  energy  to 
the  previously  non-turbulent  fluid  is  accomplished  by  the  energy 
containing  eddies,  which  transfer  energy  to  the  vorticity  diffused  into 
the  outer  fluid  by  the  small  eddies.  In  the  above  manner,  turbulence 
is  generated  into  previously  quiescent  fluid.  The  turbulent  wake 
"front"  represents  a  boundary  between  fluid  which  possesses  vortic¬ 
ity  and  vorticity  fluctuations,  and  fluid  which  does  not.  Corrsin  and 
10 

Kistler  describe  the  front  as  a  "laminar  superlayer",  of  very 
small  thickness,  across  which  vorticity  decreases  from  values 
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characteristic  of  the  turbulent  flow,  to  very  low  values.  Vorticity 
is  communicated  locally  to  new  fluid  across  this  thin  layer,  prin¬ 
cipally  by  the  smallest  eddies  of  the  turbulent  flow,  resulting  in  a 
'propagation"  of  the  front  into  new  fluid.  The  new  fluid  thus  acquires 
vorticity,  and  vorticity  fluctuations,  but  it  acquires  no  large  scale 
random  motions  by  the  above  mechanism.  Consequently,  the  newly 
engulfed  fluid  would  tend  to  remain  close  to  its  non-turbulent 
"trajector/",  under  the  effect  of  the  primary  turbulence  prop,  gation 
me  lanism.  It  would  therefore  not  be  really  mixed  as  yet,  and  the 
small  scale  motions  induced  in  it  would  not  appreciably  aid  the 
molecular  mixing  process. 

If  there  were  no  other  rr’xing  mechanism,  each  additional 
"layer"  of  engulfed  fluid  would  come  in  contact  with  only  the  preced¬ 
ing  la''^er  of  engulfed  fluid,  and  effective  mixing  would  not  occur. 
Effective  mixing  would  require  a  mechanism  which  swept  aw’ay 
newly  engulfed  fluid  and  mixed  into  the  "deeper"  parts  of  the  flow, 
while  exposing  new  layers  to  fluid  bi.^ught  up  from  those  deeper 
parts.  Such  a  mechanism  does  indeed  exist,  but  it  is  quite  distinct 
from  the  primary  entrainment  mechanism,  and  is  due  to  the  turbu¬ 
lent  motions  on  another  scale.  Indeed,  in  addition  to  the  small  scale 
eddies,  the  turbulent  flow  exhibits  larger  scale  eddy  motions,  which 
convect  turbulent  fluid  from  one  part  of  the  wake  to  another.  (The 
large  eddies  also  serve  to  contort  the  turbulent  front  and  thereby 
greatly  increase  the  rate  of  entrainment  of  new  fluid,  since  the 
latter  depends  on  the  total  area  of  the  "laminar  superlayer''.  )  The 
newly  engulfed  fluid  will  not  become  mixed  into  the  core  until  it  is 
convected  by  one  of  those  larger  eddy  motions.  The  scale  of  the 
inhomogeneities  of  density,  temperature,  etc.  introduced  into  the 
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turbulent  core  in  this  manner  then  depends  on  how  much  new  fluid  is 
engulfed  by  the  propagation  of  the  turbulent  front.  It  may  be  roughly 
estimated  as  follows:  7'he  main  convectiv'e  action  of  the  turbulence 
is  carried  on  by  the  larger  eddies  of  the  turbulent  flow,  of  scale  L, 
say.  The  time  required  by  those  eddies  to  sweep  away  nevvly 
engulfed  fluid  is  of  the  order  T  ~  (L/u'),  where  u'  is  the  mean 
velccity  of  the  larger  eddies.  Over  a  period  of  time  T  the  turbulent 
front  will  propagate  a  distance  D  ~  V  T,  where  V  is  the  mean  prr  ■ 
pagation  velocity  of  the  front.  The  "transverse"  scale  of  the  en¬ 
gulfed  fluid  will  be  of  the  order  of  D,  while  its  "longitudinal"  scale 
(along  the  boundary)  will  be  some  fraction  of  L.  These  scales  m.ay 
now  be  compared  with  the  scales  of  inhomogeneities  at  which  mole¬ 
cular  diffusion  becomes  dominant.  The  latter  scale  is  comparable 
to  the  dissipation  or  Kolmogorov  scale  when  the  kinematic  viscosity 
and  diffusivity  are  comparable.  In  a  well  developed  turbulent  flow, 
it  can  be  expected  to  be  two  or  more  orders  of  magnitude  below  the 
macroscale  of  turbulence,  namely  L.  The  question  remains  whether 
D  is  also  larger  than  the  Kolmogorov  scale.  From  the  above  dis¬ 
cussions,  we  may  write 


*  y"' 

D~  V  r~L-V 
u 


(1) 


where  Vo, is  the  free  stream  velocity,  and  AV  is  the  velocity  differ¬ 
ence  across  the  wake.  The  rtUio  (Va>/AV)  of  the  free  stream  ’-■-iocity 

2 

to  the  wake  velocity  defect  can  be  roughly  taken  u.s  (V^/ ^V)  ~  8y^ 

{pi poo)  where  y^  is  the  non-dimensional  wake  radius,  normalized  with 
respect  to  the  body  diameter  and  {pf  Pm)  is  the  ratio  of  mean  wake  to 
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free  stream  density.  The  ratio  (AV/u')  is  approximately  2.  accord¬ 
ing  to  Subsonic  turbulent  wake  results.  The  ratio  (V  is  approxi¬ 

mately  the  same  as  the  wake  "slope”  (dyi'/dz).  Assuming  an  approxi¬ 
mate  (non-dimensional)  wake  growth  law  of  the  form  ~  0.  , 

where  z  is  downstream  distance  expressed  in  body  diameters,  the 
above  estimates  can  be  substituted  in  Eq.  (1),  which  yields: 


(2) 


In  the  region  beyond  the  first  ten  or  twenty  body  diameters 
where  the  present  model  can  be  expected  to  apply,  the  mean  density 
p  in  the  wake  will  not  be  less  than  a  few  times  below  the  free  stream 
density,  so  that  D  will  be  of  the  order  of  one -tenth  of  L,  which  may 
reasonably  be  assumed  to  be  significantly  greater  than  the  Kolmogorov' 
scale.  Thus  it  appears  reasonable  to  assume  that  the  combined 
effects  of  the  propagation  of  the  turbulent  front  and  the  convective 
action  of  the  large  eddies  of  the  turbulence  tend  to  introduce  into  tlie 
wake  fluid  elements  of  a  scale  large  compared  to  the  Kolinogorov  or 
dissipative  scale  of  the  turbulence,  so  that  these  must  be  reduced  in 
scale  by  the  cascading  process  of  turbulence  before  effective  molecu¬ 
lar  mixing  can  occur.  The  relative  time  scales  for  the  turbulent 
cascade  process  and  molecular  diffusion  or  conduction  are  discussed 
in  Appendix  B.  It  is  shown  in  that  Appendix  that  for  large  Reynolds 
numbers,  the  critical  assumption  that  the  turbulent  cascade  time  is 
short  compared  tc  .he  diffusion  time  scale  is  valid. 

In  summary,  the  engulfment  of  the  outer  wake  by  the  turbulent 
core  may  be  pictured  as  a  process  whereby  relatively  large  fluid 
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elementb  (i.  e.  ,  large  compared  to  the  Kolmogorov  scale)  are  con- 
vected  into  the  turbulent  core  by  the  action  of  large  scale  eddies, 
while  at  the  same  time  they  are  distorted  by  the  shearing  effects  of 
small  scale  eddy  motions until  the  scale  of  the  inhomogeneities 
has  been  reduced  (and  the  gradients  steepened)  to  the  point  where 
molecular  diffusion  erases  those  inhomogeneities  and  mixes  the 

vO 

engulfed  fluid  with  the  core  fluid  on  a  molecular  scale. 

As  a  result  of  the  mixing  process  described  above,  the  wake 
can  be  expected  to  consist  essentially  of  two  distinguishable  com¬ 
ponents,  namely  a  component  which  has  not  been  mixed  a  mole¬ 
cular  scale  with  the  ’’main"  flow,  and  whose  properties  are  close  to 
the  properties  it  would  have  if  it  were  not  engulfed  in  the  core,  and 
a  more  or  less  homogen.'^ous  component  which  consists  of  the  initial 
core  fluid  augmented  by  that  portion  of  the  engulfed  fluid  which  has 
been  molecularly  mixed.  Since  both  the  entrainment  process  and  the 
mixing  process  are  random  in  nature,  the  above  description  of  the 
wake  must  be  assumed  to  apply  'in  the  mean,  '  so  to  speak. 

The  model  for  the  wake  structure  formulated  below  assumes 
that  there  is  no  significant  fluid  component  v  hich  is  neither  unmixed 
nor  corrnletely  (miolecularly)  mixed.  This  assumption  is  more 
important  when  chemical  reactions  are  concerned,  in  which  case  the 
times  characteristic  of  chemical  reactions  must  be  long  compared 


ULi 

n* 

In  the  absence  of  molecular  diffusion  and  conduction,  the  effect  of 
eddy  diffusivity  (or  random  convection)  is  to  break  up  large  volume 
elements  into  smaller  ones  (down  to  the  cut-off  scale)  and,  there¬ 
fore,  to  increase  mean  gradients  rather  than  decrease  them.  The 
incre  ise  in  mean  gradients  leads,  however,  to  increased  molecular 
dissipation,  so  that  eddy  diffusion  indirectly  leads  tu  dissipation  of 
inhomogeneities . 
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to  the  time  required  for  fluid  elements  to  jump'  across  the  Kolmo¬ 
gorov  scale.  However,  the  basic  concept  used  in  the  model,  namely 
that  a  macroscopic  engulfment  and  mixing  precedes  molecular  mix¬ 
ing,  is  not  dependent  upon  that  assumption. 

In  the  model  of  wake  structure  described  above,  the  spatial 
variations  in  wake  properties,  such  as  mass  density,  electron  den¬ 
sity,  and  temperature,  are  regarded  as  arising  from  the  existence 
of  'unmixed'  fluid  elements  intermingled  with  'mixed'  fluid,  as  well 
as  from  inhomogeneities  which  may  exist  within  or  between  unmixed 
fluid  elements.  Inhomogeneities  in  '^^he  mixed  component  are  assumed 
small,  and  will  in  fact  be  neglected  in  the  following. 

The  mean  density  fluctuations  which  can  be  predicted  on  the 

basis  of  a  two-component  wake  model  such  as  the  one  described 

above  are  in  semi -quantitative  agreement  as  will  be  shown  below, 

^  *7  12 

with  the  large  mass  density  fluctuations  apparently  found  ’  ’ 
in  schlieren  measurements  of  density  fluctuations  in  the  wakes  of 
hypervelocity  pellets  at  large  downstream  distances  from  the  pellets. 
The  apparent  magnitude  of  the  measured  fluctuations  at  large  down¬ 
stream  distances,  of  the  order  of  at  least  one  tenth  at  an  axial  dis¬ 
tance  of  several  hundred  body  diameters,  would  seem  to  rule  out  the 
possibility  that  the  fluctuations  are  created  within  the  turbulent  flow 
by  compressibility  and  dissipative  effects.  Indeed  it  is  shown  in 
Appendix  A  that  such  effects  would  lead  to  considerably  smaller 
fluctuations.  Similarly,  interaction  between  the  turbulence  and  the 
wake  chemistry  must  be  ruled  out  for  the  relatively  low  tempera¬ 
tures  in  the  observed  ballistic  range  wakes. 
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In  the  far  wake,  where  the  turbulent  velocities  have  presum¬ 
ably  died  down,  the  above  model  is  not  expected  to  hold.  The  eddy 
diffusivity  is  probably  no  longer  very  important,  and  the  break  up  of 
the  lumps  no  longer  occurs  very  efficiently.  The  temperature 
inhomogeneities  then  decay  essentially  by  conduction.  In  fact,  suf¬ 
ficiently  far  downstream  all  velocities  will  have  essentially  died 
down.  Thus,  the  final  steps  of  thermal  equalization  proceed  purely 
by  molecular  diffusion.  The  limit  of  application  of  the  two-fluid' 
model  is  undoubtedly  dependent  on  Reynolds  number  and  body  geom¬ 
etry,  and  it  may  not  be  applicable  at  all  when  the  turbulent  intensity 
is  very  low,  corresponding  to  small  Reynolds  numbers. 

The  model  may  also  require  modification  for  application  to 
the  wakes  of  a  e-entry  bodies,  particularly  close  to  the  body,  as  the 
character  of  the  turbulence  in  the  near  wake  may  be  quite  different, 
and  compressibility  effects  in  particular  may  be  important. 
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3.  WAKE  MIXING  EQUATIONS 


A  simple  set  of  equations  describing  a  quasi  one -dimensional 
wake  with  a  mean  turbulent  boundary,  and  incorporating  the  physical 
features  of  wake  mixing  discussed  in  the  previous  section,  is  pre¬ 
sented  below.  The  equations  are  similar  to  those  presented  in 
References  4  and  5.  As  in  those  references,  the  growth  of  the 
(mean)  turbulent  wake  is  represented  by  the  flux  of  inviscid  fluid 
into  the  turbulent  core  through  the  wake  boundaj  /.  The  entering 
fluid  then  forms  part  of  the  unmixed  portion  of  the  core  for  a  mean 
distance  ^  (measured  in  body  diameters),  after  which  it  is  homoge¬ 
neously  mixed  with  the  remainder  of  the  wake  by  molecular  effects. 
The  actual  molecular  mixing  of  fluid  entering  at  a  particular  station 
in  the  wake  need  not  occur  suddenly  after  a  distance  but  may  be 
spread  over  a  (small)  range  of  lags  centered  about  4*  This  is  dis¬ 
cussed  further  in  the  following  section. 

The  turbulent  wake  is  then  characterized  by  the  pressure 
(assumed  constant  across  the  wake),  and  by  the  velocity,  density, 
enthalpy  and  chemical  composition  (if  chemistry  is  included)  of  the 
homogeneously  mixed  and  xtnmixed  portions  of  the  wake,  together 
with  the  relative  fractions  of  each,  as  a  function  of  downstream 
distance.  The  values  of  the  velocity,  density,  enthalpy,  and  chemi¬ 
cal  composition  of  the  inviscid  gas  at  the  edge  of  the  turbulent  core 
must  be  assumed  to  be  specified.  The  wake  evolution  can  then  be 
described  by  the  one -dimensional  equations  of  conservation  of  mass, 
momentum  and  energy,  by  an  equation  describing  the  mixing  lag,  and 
by  the  equations  expressing  the  chemical  reactions  or  the  thermo¬ 
dynamic  characteristics  of  the  gas  (equations  of  state)  if  chemical 
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reactions  can  be  ignored.  In  the  latter  case,  v/hich  is  applicable  to 
wakes  of  relatively  low- speed  projectiles  in  ballistic  ranges  for 
instance,  the  conservation  equations  specifying  the  quasi  one- 
dimensional  model  can  be  expressed  as  follows  (z  denotes  down¬ 
stream  distance): 

Conservation  of  Mass: 


dz 


(p  V  dA)  +  3—  (p  V  A  ) 
'^u  u  '  dz  ^m  m  m 


=  P-V.  — 
*^1  i  dz 


(3) 


Conservation  of  Momentuir; 


dz 


(p  dA)  +  |-  (p  A  )  - 
^u  u  dz  m  m  m 


Z  dA 
P  V.  3— 
i  dz 


(4) 


Conservation  of  Energy: 


P  V 

m  m 


+ 


=  p.V. 

i  i  dz 


(5) 


The  foregoing  equations  are  the  usual  conservation  equations  for 
one imensional  flow  in  a  channel  with  mass  ■'.ddition  at  the  boundary 
(at  the  same  pressure  as  in  the  channel),  for  a  flow  consisting  of 
two  components.  In  the  foregoing  equations,  p^  is  the  density  of  the 
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inviscid  gets  at  the  edge  of  the  turbulent  core,  V.  is  its  velocity,  H 

2  -  ^ 

its  total  enthalpy,  H  =  [h.  +  (V^/2)],  and  p  is  the  pressure.  The 
density  p,  velocity  V,  and  static  enthalpy  h,  for  the  mixed  compo¬ 
nent  and  for  the  unmixed  portions  of  the  turbulent  wake  are  denoted 

by  the  subscripts  m  and  u  respectively.  The  areas  A  and  A  are 

m  u 

the  'partial  areas'  of  hot  and  cold  fluid,  and  their  sum  is  equal  to 
the  wake  area  A: 


A  +  A  =  A 
m  u 


(6) 


In  the  above  cited  references,  it  had  been  assumed  that  the 
unmixed  fluid  could  be  characterized  by  a  single  average  value  of  p 
and  V.  This  assumption  is  not  necessary  in  the  calculations  pre¬ 
sented  below,  in  which  the  properties  of  the  unmixed  fluid  component 
are  assumed  known. 

The  mixing  lag  can  be  introduced  in  the  mixing  model  by  an 
equation  of  the  form 

A  )  =  p.(z  -  OV,(z  -  (7) 

dz  m  m  m  i  i  dz 

where  ^  is  the  lag  distance  defined  above,  which  expresses  the 
assumption  that  fluid  ingested  at  station  (z  -  4)  becomes  part  of  the 
mixed  component  at  station  z,  that  is,  with  a  lag 

The  solution  to  the  wake  mixing  equations,  Eqs.  (3)  to  (7) 
becomes  particularly  simple  if  it  is  assumed  that  the  unmixed  fluid 
component  is  entirely  unaffected  in  its  thermodynamic  and  chemical 
properties  by  the  random  convection  to  which  it  is  subjected  within 
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the  turbulent  core.  Only  the  mixed  fluid  properties  then  remain  to 
be  specified.  Equation  (7)  then  shows  that  in  computing  the  proper¬ 
ties  of  that  component,  the  wake  mixing  model  is  equivalent  to  an 

13 

instantaneous  mixing  model  ,  but  with  an  effective  mixing  boundary. 

The  relation  of  the  effective  mixing  boundary  to  the  mean  turbulent 

boundary  is  depicted  in  Fig.(l).  Inviscid  fluid  entering  the  turbulent 

boundary  at  (zj.,  yj.)  survives  molecular  mixing  for  a  distance  i,  so 

that  it  would  have  reached  the  point  (z  ,  y  )  if  it  had  continued 

m  m 

along  its  laminar  trajectory.  Since  the  random  convection  by  the 

turbulent  velocity  field  in  the  wake  is  assumed  not  to  affect  the  fluid, 

it  is  effectively  mixed  at  (z  ,y  ).  Thus,  the  mixing  boundary  can 

be  constructed  from  the  mean  wake  boundary  by  simply  following 

fluid  elements  on  an  extension  of  the  laminar  streamlines  from  their 

point  of  crossing  of  the  turbulent  boundary  to  a  point  ^  downstream. 

The  mount  of  unmixed  fluid  is  then  that  contained  between  y  and  y  . 

f  m 

The  wake  mixing  equations  describing  the  mixed  fluid  com¬ 
ponent  are  then  simply  the  conservation  equations  for  one -dimensional 
channel  flow,  with  or  without  chemical  reactions,  and  with  mass 
injection  at  the  effective  mixing  bound^-ry.  In  order  to  solve  those 
equations,  the  lag  must  be  specified,  and  it  is  discussed  in  the 
following  section. 
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f9-Zl  ‘BHH 


I'l 


TURBULENT  FRONT 


4  MIXING  LAG  DISTANCE 

In  Eq.  (7),  the  asi^amed  lag  in  the  molecular  mixing  of 
ingested  fluid  into  the  core  is  represented  by  a  lag  distance  4-  This 
distance  is  approximately  equal  to  some  lag  time  T  multiplied  by  the 
mean  wake  velocity.  The  lag  time  T  is  clearly  related  to  the  time 
scale  for  dissipation  of  inhomogeneities  in  a  turbulent  flow.  This 
lag  is  determined  by  the  rate  of  molecular  mixing  in  a  turbulent 
flow  field,  and  is  therefore  related  to  the  rate  of  destruction  of 
scalar  inhomogeneities  in  such  a  flow  field.  The  only  such  rate 
presently  available  from  the  theory  of  turbulence  is  that  for  the 
mixing  of  a  passive  scalar  in  a  homogeneous  steady  state  field  of 
turbulence  ’  and  even  then  the  rate  is  only  known  at  large  values 
of  the  time,  when  the  effect  of  the  initial  distribution  of  scalar 
inhomogeneities  has  disappeared.  In  that  case,  the  (logarithmic) 
rate  for  the  dissipation  of  scalar  inhomogeneities  is  the  same  as  the 
dissipation  rate  of  kinetic  energy.  There  is,  of  course,  no  good 
a  priori  reason  to  believe  that  the  distribution  of  inhomogeneities  in 
the  wake  mixing  problem  would  be  the  same  as  that  in  such  a  homo¬ 
geneous  steady-state  field  of  turbulence.  However,  it  is  not  unrea¬ 
sonable  to  hope  that  the  mean  rate  of  dissipation  of  inhomogeneities 
predicted  for  steady- state  homogeneous  turbulence  will  not  be  too 
far  off  in  the  present  problem. 

The  time  scale  for  dissipation  of  inhomogeneities  in  homo¬ 
geneous  turbulence  can  be  estimated  by  using  a  relation  due  to 
Onsager^^  describing  the  process  of  turbulent  cascade  of  energy 
from  large  scale  to  small  scale  eddies: 
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(8) 


dk  dt 


where  a  €  k  is  the  ‘'lifetime"  of  an  eddy  of  wave  number  k. 
Eq.  (8)  describes  an  accelerated  transfer  of  energy  from  large  to 
small  scale  eddies.  In  that  equation,  c  is  the  rate  of  dissipation  of 
kinetic  energy  of  turbulence  by  viscous  effects. 

Equation  (8)  can  be  integrated  ‘o  o'  ain  the  characteristic 
dissipation  time  scaler  required  for  the  breakup  of  an  inhomogeneity 
of  characteristic  scale  D  =  down  to  the  dissipation  or  Kolmo¬ 
gorov  scale  ri  -  assumption  that  the  dissipation  of 

scalar  inhomogeneities  can  be  described  in  terms  of  the  turbulent 
cascade  process.  Thus,  Eq.  (8)  yields 


T 


T-T 


(9) 


Since  the  Kolmogorov  scale  k  is  much  larger  than  k_,  the 

j  3 

term  k  can  be  neglected  in  the  above  equation.  P^urthermore , 

Q 

the  dissipation  rate  €  can  be  written  approximately  as: 


3  du'^  _  3  u'^ 
2  dt  “  2  L 


(10) 


where  3/2  u'  is  the  mean  turbulent  kinetic  energy,  and  L  is  the 
scale  of  the  energy  containing  eddies  of  turbulence.  Using  the  above 
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equation,  and  assuming  that  the  initial  scale  of  the  inbomogeneities 
is  given  by  Eq,  {2.),  Eq.  (  9  )  can  oe  rewritten,  approximately; 


I  du'^ 

,2  dt 
u 


dt 


(11) 


where  p  is  a  mean  wake  velocity  at  the  initial  time. 

Substituting  downstream  distance  for  time,  Eq.  (11)  yields; 


1  du'^ 

dz 

u 


dz 


(12) 


where  p  is  to  be  evaluated  at  (z  -  Q. 

Finally,  the  turbulent  velocity  can  be  expressed  in  terms 
of  the  mean  wake  and  full  stream  velocities,  by  writing 


(13) 


Let 


fi(z) 
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and 


=  G(z) 


Then  Eq.  (12)  can  be  immediately  integrated  to  give: 


G(z)  =  G(z  - 


^  n(z) 


(14) 


The  function  n(z),  which  is  the  ratio  of  turbulent  to  mean  flow  kinetic 
energy  m  the  wake,  tends  in  subsonic  wakes  to  a  constant  value  at 
large  downstream  distances,  of  the  order  of  a  few  hundred  body 
diameters  for  two-dimensional  wakes  and  probably  much  sooner  for 
axi-symmetric  wakes.  On  the  o..her  hand,  at  small  downstream 
distances,  at  which  Q(z)  decays  fairly  rapidly,  ^  is  a  small  quantity, 
compared  to  the  decay  scale  of  n(z).  Thus,  the  ratio  [P'  •  -  ^)3/[n(z)] 
will  be  approximately  constant  and  equal  to  unity  everywhere  in  the 
wake.  With  this  assumption,  Eq.  (14)  becomes  a  simple  implicit 
equation  for  which  depends  essentially  only  on  the  form  of  the 
function  G(z)  representing  the  decay  of  mean  wake  velocity  wuth 
downstream  distance.  It  is  not  sensitive  to  the  exact  value  of  either 
the  macroscale  or  the  dissipaticn  scale  of  the  turbulence,  nor  to  the 
exact  form  of  the  function  n(z). 

If  a  DOwer  law  decay  is  assumed  for  the  w'ake  velocity,  of  the 
form  G  (z)  =  const,  x  z  ,  Eq.  (14)  yields 
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(15) 


This  form  has  been  fitted  in  the  numerical  calculations  presented 
below . 

Clearly  one  cannot  expect  Eq.  (14)  or  (15)  to  provide  an 
exact  estimate  of  the  lag.  In  fact,  .uch  an  exact  estimate  is  not 
very  meaningful,  since  ^  is  intended  to  represent  some  mean  mixing 
lag,  with  the  actual  mixing  of  any  giv-^en  fluid  element  occurring  over 
some  range  of  downstream  distances.  Such  a  range  in  mixing  dis¬ 
tance  couid  be  expected  to  arise  from  the  fact  that  the  survival  time 
T  of  a  fluid  element  will  exhibit  a  range  of  values,  from  the  fact  that 
a  given  element  will  not  mix  entirely  at  one  time,  and  finally  from 
the  fact  that  there  is  a  velocity  profile  in  the  w'ake,  leading  to  a 
range  in  ^  even  for  a  fixed  survival  •’ime.  However,  if  that  range  is 
fairly  small  compared  to  the  length  scale  over  which  the  outer 
(inviscid)  field  properties  vary  appreciably,  a  "sharp”  lag  ^  can  be 
used  to  represent  mixing  of  a  fluid  element  over  a  range  of  distance 
after  it  is  ingested.  Indeed,  a  more  general  form  of  Eq.  (7)  rep¬ 
resenting  mixing  over  a  range  of  downstream  distance  might  be 
given  as 


The  use  of  a  lag  equation  summing  the  contributions  to  the  mixed 
flow  at  a  given  station  from  the  unmixed  flux  at  all  preceding 
stations  was  suggested  by  Dr.  A.  Nayfeh  of  Heliodyne,  though  in  a 
somewhat  different  form. 
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z 


jL 

dz 


V  A  \  =  f 

m  J 


z 


0 


pJz')V.{z') 


d  A(z ' )  50  (z,  z ' ) 
dz'  5z 


(16) 


where  0{z,z')  is  the  portion  of  the  inviscid  flux  per  unit  area  enter¬ 
ing  the  core  at  z'  which  has  survived  to  station  z,  and  therefore 
50/ dz  is  a  rate  of  mixing  at  station  z  of  inviscid  fluid  which  entered 
the  core  at  z'.  (Obviously  0  =  0  when  z  <  z'.  ) 

The  discussion  of  mixing  in  a  turbulent  flow  presented  in  the 
preceding  sections  suggests  that  50/ 5z  is  a  functio  ,  which  peaks  at 
some  particular  value  of  z  -  z',  namely  when  the  shredding  effect  of 
turbulence  has  reduced  the  scale  of  the  inhomogeneiti^s  to  the  Kol¬ 
mogorov  scale,  and  is  small  everywhere  else. 

Thus,  particularly  in  the  case  of  slender  body  wakes  (where 
the  inviscid  flow  field  does  not  vary  much),  but  without  much  error 
in  all  cases  where  chemical  reactions  do  not  make  the  flow  proper¬ 
ties  vary  too  rapidly  with  distance,  we  may  rewrite  Eq.  (16)  as 


|-/p  V  A  )  ^  p, 
dz  y  m  m  mi  l 


;(Z  -  0V.(Z 


dA  (z  -  i)  f  50  ,  , 

-d(.  -5,  ■  J  sT 

^0 


,  dA  (z  -  6) 

P.(z  -  4)V,(z  -  C)  — 


(17) 
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where  |  is  now  the  value  about  which  (d0/Sz)  peaks.  If  the  peak  of 
(30/92)  Is  assumed  infinitely  sharp, 


30 

-  4) 


6(z  -  i)  . 


and 


3z 


d(z  -  €) 

dz 


6(z  -  i)  , 


confirming  Eq.  (7).  Thus  the  wake  mixing  formulation  using  a 
sharp  lag  can  be  used  to  describe  gradual  mixing  over  a  range  small 
compared  to  the  range  over  which  the  inviscid  flow  properties 
change  significantly. 
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5.  NUMERICAL  RESUI.TS  AND  DISCUSSION 


In  the  present  section,  a  numerical  solution  of  the  wake 
mixing  equations  using  a  simplified  set  of  an  effective  mixing  bound- 
ar/  is  presented  and  discussed.  It  should  be  stressed  at  the  outset 
that  the  model  is  quite  crude,  and  that  in  particular  it  does  not 
account  very  well  for  what  might  be  termed  dynamic  effects,  such 
as  drag  between  the  mixed  and  unmixed  components,  the  effects  on 
pressure  of  mixing  of  the  two  components,  and  so  forth.  It  is  pri¬ 
marily  what  its  name  implies,  i.  e.  ,  a  mixing  model,  and  is  intended 
for  application  to  reacting  wakes.  However,  the  low  speed  non- 
reaciing  case  is  of  interest,  particularly  in  comparing  its  predic¬ 
tions  with  experimental  results.  The  calculations  were  made  to 
correspond  to  the  case  of  a  one  centimeter  diameter  body  traveling 
at  some  10,000  feet/sec. 

The  mean  wa  e  boundary  used  in  the  numerical  computations 
was  obtained  from  a  -j  power  law  fitted  to  experiment.  The  lag 
function  was  approximated  by  a  linear  dependence  of  ^  on  down¬ 
stream  distance  over  several  wake  segments,  with  the  coefficients 
chosen  to  roughly  match  Eq.  (15).  The  value  of  n  in  Eq.  (15)  over 
the  various  segments  was  estimated  by  reference  to  measured  data 
on  velocity  decay  in  wakes and  {pi p^)  was  estimated  from  a 
theoretical  calculation  by  L.  Hromas^^.  From  such  estimates,,  the 
lag  V  as  assumed  to  rise  from  0.  05z  for  z  less  than  50  body  diam¬ 
eters  to  0.  Iz  between  30  and  150  diameters,  and  increasing  to  0.  5z 
for  z  greater  than  300  diameters,  with  additive  constants  used  to 
make  ^  continuous.  Beyond  500  diameters,  ^  was  approximated  m 
reach  a  constant  value  corresponding  to  a  lag  time  of  approximately 
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two  milliseconds,  which  is  approximately  the  diffusion  time  for 
molecular  diffusion  from  a  one  centimeter  diameter  sphere  at  a 
temperature  of  approximately  SOO'^K  and  six  percent  atmospheric 
pressure. 

The  specification  of  the  inviscid  flow  outside  the  turbulent 

core  was  obtained  from  a  hemisphere  cylinder  solution  also  supplied 

by  L-  Hrom.as^^.  The  free  stream  velocity  in  the  numerical  case 

used  was  about  9600  ft/  s  and  the  free  stream  density  used  was 
_  2 

(p^/p^)  ~  6  X  10  corresponding  to  an  attitude  of  bO,  000  ft.  An 
equilibrium  gas  was  assumed  in  the  numerical  computations. 

The  results  obtained  by  solving  the  one-dimensional  channel 
flow  with  an  effective  mixing  boundary  are  shown  in  Figs.  (2)  to  (4). 

Figure  (2)  shows  the  variation  of  the  density  p  along  the 

wake  as  welt  as  the  mean  density  p^  of  the  unmixed  component.  It 

can  be  seen  that  p  becomes  subs'^antially  greater  than  p  as  we 

u  m 

proceed  downstream.  This  leads  to  a  relatively  large  ratio 
and  the  possibility  of  large  relative  fluctuations,  as  discussed  in  sub¬ 
section  2-4.  The  ratio  T  is  plotted  in  Fig.  (3)  along  with  the  unmixed 

P 

volume  fraction  c. 

Finally,  the  relative  root-rrean  square  gas  density  fluctuation 
{p' / p)  is  plotted  in  Fig.  (4).  The  gas  density  fluctuation  is  obtained 
from  the  equation 


jE-  =  J- 

P  P 


(18) 


which  accounts  for  fluctuations  "xising  from  inhomogeneity  of  the 
unmixed  gas. 
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DOWNSTREAM  DISTANCE,  X/D 

I  of  ratio  of  mixed  fluid  density  to  mean  unmix 
mixed  volume  fraction  e  with  downstream  dist 
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Fig.  4  Normalized  mass  density  fluctuations  along  wake  according  to  mixing 
with  lag  model.  The  experimental  values  shown  have  been  recently 
revised  sharply  downward  by  the  experinjenters  12,  although  precise 
new  values  have  not  yet  been  made  available. 


The  experimental  results  obtained  at  Lincoln  Laboratory 
using  schlieren  methods  agree  with  tliC  predicted  trend  in  the  data, 
with  a  maximum  observed  amplitude  around  0.  I,  for  pellets  travel¬ 
ing  at  approximately  8500  ft/sec.  There  appears  to  be  some  un¬ 
certainty,  however,  in  the  determination  of  the  absolute  magnitude 
of  the  density  fluctuations  which  they  measure. 

The  present  paper  has  been  based  upon  a  characterization  of 
wake  mixing  in  terms  of  propagation  of  turbulence  by  small  scale 
eddies,  random  convection  of  fluid  elements  wuthin  the  turbulent 
core  by  large  scale  eddies,  and  dissipation  of  inhomogeneities  by 
the  combined  action  of  distortion  of  fluid  elements  by  turbulent 
motions  and  molecular  diffusion.  Such  a  characterization  is  consis¬ 
tent  with  the  theory  of  (low-speed)  turbulence  as  it  exists  today.  It 
leads  to  a  two-component  wake  model,  in  which  one  component  con¬ 
sists  of  fluid  which  is  relatively  homogeneous  and  well  mixed,  while 
the  other  consists  of  fluid  newly  entrained  into  the  turbulent  core, 
and  therefore  retaining  its  unmixed  properties.  In  the  above  model, 
the  time  delay  in  thorough  (molecular)  mixing  of  the  latter  compo¬ 
nent  was  represented  by  a  lag  time  or  lag  distance. 

A  set  of  one -dimensional  equations  derived  for  such  a  model 
was  developed,  and  the  solution  to  a  simplified  version  of  those 
equations  for  a  relatively  low'  speed  v/ake  provided  predictions  of 
density  fluctuations  in  the  wake  in  good  agreement  with  available 
experimental  data. 

The  model  developed  above  is  intended  for  use  in  predicting 
the  properties  of  turbulent  wakes  of  high  speed  re-entry  vehicles. 
For  such  wakes,  the  chemical  reactions  in  the  w'ake  are  of  primary 
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importance,  and  the  model  requires  extension  and  further  refinement 
in  order  to  provide  reliable  predictions  for  reacting  wakes  of  high 
Speed  bodies,  though  it  can  be  used  for  rough  calculati  )ns  by  includ¬ 
ing  reactions  in  the  conservation  equations  for  each  specie  of 
importance.  In  any  event,  the  ability  of  the  model  to  predict  some 
of  the  features  of  wakes  of  low-speed  bodies  is  an  encouraging  indi¬ 
cation  of  the  ''^alidity  of  the  concept  of  a  two-component  wake. 

Finally,  the  validity  of  the  model  could  be  more  critically 
tested  if  measurements  of  both  gas  and  electron  density  fluctuations 
in  the  wake  were  available,  as  the  model  makes  definite  predictions 
concerning  the  relation  between  those  fluctuations. 
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APPENDIX  A 


In  this  appendix,  a  rough  estimate  of  the  magnitude  of  density 
oscillations  which  can  be  generated  in  a  turbulent  wiike  by  dissipation 
and  compressibility  effects  is  provided. 

The  order  of  magnitude  of  the  density  fluctuations  that  could 
be  generated  in  a  turbulent  wake  may  be  estimated  as  follows: 
assume  that  the  source  of  the  fluctuations  j.s  dissipation  or  pV  work 
and  is  generated  by  the  fluctuating  velocities  rather  than  convection 
of  cold  high  density  gas  into  hot  low  density  gas.  Since  there  is  no 
way  to  support  large  pressure  fluctuations  in  a  free  flow  like  a  w'ake, 
beyond  a  few  downstream  diameters  (particularly  for  free  stream 

4 

velocities  of  10  ft/  sec  or  so),  we  assume  that  the  fluctuations  occur 
at  constant  pressure.  For  the  case  of  fluctuations  generated  by 
viscous  dissipation,  it  may  be  assumed  that  at  most  all  the  random 

kinetic  energy  per  unit  mass  in  a  given  unit  mass  of  fluid  will  be 

'I' 

converted  into  heat  in  generating  a  density  change. 

The.  kinetic  energy  E,  ,  per  unit  mass  in  the  random 
°  turb 

velocity  field  is  approximately 


where  u'  is  the  root  mean  square  turbulent  velocity,  AV  is  the  wake 
velocity  defect,  and  V^is  the  free  stream  velocity. 


This  is  an  upper  limit  by  any 


standard. 
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The  ratio  (u*  /  AV  )  can  be  e5timated  to  be  roughly  one-tenth 
bey  ond  a  lew  tens  of  body  die  meter.  The  approximate  relation  (see 
Section  2.  1) 


(A-2) 


may  be  used  for  the  wake  v  ocity  decay. 


On  the  other  hand,  assuming  that  all  the  kinetic  energy  of 
nee  is  con' 
density  change  is 


turbulence  is  converted  to  heat,  so  that  c  AT  ~  E  ,  ,  the  resulting 

P  turb  ^ 


E 

AP  _  ^  _  turb 

p  ‘  T  ■  ~  c  T 
P 


(A-3) 


where  h  is  the  local  static  enthalpy. 

2 , 

The  ratio  (V^/h)  can  be  estimated  by  writing  it  in  the  form 


h 


^4-4) 


Finally,  substituting  Eqs.  (A-1),  (A-2)  and  tA-4)  in^^o  (A-3),  and 

1/3 

assuming  ~  (l/2)z  ,  the  density  fluctuation  can  be  estimated  as 


Ia£|  ^ 

p 


(A-5) 
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where  p  is  the  mean  wake  density.  From  the  above  equation,  it  is 
clear  that  at  downstream  distances  of  the  order  of  one  to  several 
hundred  body  diameters,  the  density  fluctuations  due  to  dissipation 
can  be  expected  to  be  of  the  order  of  1%  or  less.  The  same  esti¬ 
mates  can  be  shown  to  result  if  it  is  assumed  that  the  density  fluctua¬ 
tions  are  generated  by  conversion  of  turbulence  energy  into  pV  work. 

Thus,  it  IS  safe  to  assume  that  if  fluctuations  in  density  of 
the  order  of  several  tens  percent  do  exist  in  the  wake,  as  some 
evidence  indicates,  they  must  arise  from  a  random  mixing  effect. 
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APPENDIX  B 


In  the  present  appendix  characteristic  time  scales  for  turbu¬ 
lent  cascade  or  dissipation  processes  and  for  molecular  diffusion 
are  estimated  and  the  validity  of  the  mixing  model  used  is  discussed 
in  terms  of  those  time  scales. 

The  present  mixing  model  assumes  that  the  rate  of  molecu¬ 
lar  mixing  of  fluid  ingested  into  the  turbulent  core  with  the  gases 
present  therein  is  limited  by  the  r  *te  at  w’hich  the  turbulent  velocity 
field  breaks  up  fluid  elements  from  a  relatively  large  scale  D  to  a 
scale  at  which  molecular  diffusion  is  highly  effectiv'e.  This  implies 
that  the  time  scale  for  this  turbulent  shredding  is  short  compared 
to  the  time  for  molecular  diffusion  to  erase  inhomogeneities  of 
scale  D.  Thus,  a  minimum  requirement  for  the  validity  of  the 
present  model  is: 


(B.  1) 


w'here  Tis  the  turbulence  cascade  time  andT,.^,.  is  the  molecular 

diff 

diffusion  time  scale.  The  latter  may  be  estimated  as 


T  (1)  ~  — 
diff  '  '  K  ’ 


(B.2) 


where  K  is  a  diffusion  coefficient.  If  we  choose  the  coefficient  of 
ambipolar  diffusion  in  air  for  K,  we  have; 
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{B.3) 


The  scale  D  may  be  assup'od  to  be  roughly  a  body  diameter 
for  the  purposes  of  the  present  estimate.  Substituting  the  above 
estimates  in  {B.  1),  we  can  express  it  approximately  as 


T  «  lOOD 


-1/2 


(B.4) 


The  time  T  has  been  assumed  in  the  present  paper  to  be 
roughly  the  same  as  the  turbulent  cascade  time.  The  latter  is 
approximately  L/u',  where  L  is  the  scale  of  the  energy  containing 
eddies  of  turbulence  and  u'  is  the  root  mean  square  turbulent  velocity. 
Assuming  approximately  L  to  be  equal  to  D,  and  u'  to  be  one  half  the 
mean  wake  velocity  AV,  (B.  4)  becomes; 


2 


1/2 


(B.5) 


Assuming  now  that  the  density  p  can  be  estimated  in  terms  of 
the  free  stream  density  as 
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and  using  the  same  approximate  velocity  decay  law  as  was  used  in 
Section  2  of  the  text,  namely 


we  can  express  the  inequality  (B.  5)  in  terms  of  a  free  stream 
Reynolds  number; 


R 


e 

OD 


1.1^  -  free  stream  viscosity. 


After  simple  substitutions,  the  result  is: 


R 


e 

as 


»  4Z 


(B.6) 


where  T  has  been  taken  throughout  as  some  mean  wake  temperature. 

The  above  inequality  is  easily  satisfied  at  least  to  downstream 
dist.'  nces  of  a  thousand  body  diameters  for  Reynolds  numbers  for 
which  the  wake  becomes  turbulent  close  to  th  body,  i.  e.  ,  for 
Reynolds  numbers  of  fifty  thousand  or  more. 

The  requirement  for  the  validity  of  the  mixing  model  in  terms 
of  the  relative  scales  for  turbulent  and  laminar  dissipation  is  some¬ 
what  more  stringent  than  indicated  by  Eq.  (B.  6).  Indeed,  the  turbu¬ 
lent  dissipation  time  scale  should  be  small  compared  to  the  molecu¬ 
lar  dissipation  ti  ne  scale  over  the  major  portion  of  the  dissipation 
’’history”,  so  to  speak,  of  a  fluid  element.  In  other  words  an 
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based  on 


inequality  of  the  form  (1)  must  hold  for  times  T  and  T,.,, 

diff 

scales  Smaller  than  D,  that  is,  as  the  cascade  process  is  occurring. 
It  can  be  shown  that  this  requirement  leads  to  a  required  Reynolds 
number  about  one  order  of  magnitude  greater  than  indicated  by  (B.  6). 
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